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Abstract 

Dirac showed that the existence of one magnetic pole in the universe could offer an 
explanation for the discrete nature of the electric charge. Magnetic poles appear naturally 
in most Grand Unified Theories. Their discovery would be of greatest importance for 
particle physics and cosmology. The intense experimental search carried thus far has not 
met with success. Moreover, if the monopoles are very massive their production is outside 
the range of present day facilities. A way out of this impasse would be if the monopoles 
bind to form monopolium, a monopole- antimonopole bound state, which is so strongly 
bound, that it has a relatively small mass. Under these circumstances it could be produced 
with present day facilities and the existence of monopoles could be indirectly proven. We 
study the feasibility of detecting monopolium in present and future accelerators. 
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1 Introduction 



The theoretical justification for the existence of classical magnetic poles, hereafter called 
monopoles, is that they add symmetry to Maxwell's equations and explain charge quanti- 
zation [TJ [2] . Dirac formulated his theory of monopoles considering them basically point 
like particles and quantum mechanical consistency conditions lead to the so called Dirac 
Quantization Condition (DQC), 



where e is the electron charge, g the monopole magnetic charge and we will use natural 
units h = c = 1. In this theory the monopole mass, m, is a parameter, limited only by 
classical reasonings to be m > 2 GeV [3]. In non-Abelian gauge theories monopoles arise 
as topologically stable solutions through spontaneous breaking via the Kibble mechanism 
(3]. They are allowed by most Grand Unified Theory (GUT) models, have finite size 
and come out extremely massive m > 10 16 GeV. There are also models based on other 
mechanisms with masses between those two extremes [3j [6]. All these monopoles 
satisfy Dirac's quantization condition as a consequence of the monopolar structure and 
the semiclassical quantization [TJ. 

All the attempts to discover or produce monopoles have met with failure. [3j El [7J 
IU [9j [10]. This lack of experimental confirmation has led many physicist to abandon the 
hope in their existence. A way out of this impasse is the old idea of Dirac [U [TTJ, namely, 
monopoles are not seen freely because they are confined by their strong magnetic forces 
forming a bound state called monopolium [TJl [TS] . 

Several cosmological scenarios compatible with all cosmo logical requirements [141 [151 
IT6l [T7] have been proposed [121 [HI [El [20] with the aim of making the existence of mo- 
nopolium compatible with the observation of ultrahigh-energy cosmic rays (UHECRS) 
[2"lT [22] • From these cosmological scenarios we have learned that the study of the mo- 
nopolium annihilation process provides us with information regarding the existence of 
monopoles, even if they are difficult to detect or to produce in a free asymptotic state. 
This phenomenon is not a novel feature of physics. Quark-gluon confinement describes 
the strong limit of Quantum Chromodynamics, the theory of the hadronic interactions, 
and their existence is proven by the detection of jets, showers of conventional hadrons. 
There is however a main difference between the two scenarios. In the monopolium case, 
the elementary constituents may be separated asymptotically, when they are orbiting far 
from each other, if the energy provided to the system is high enough, while in the quark- 
gluon case this is not possible. In practice, however, there is no big difference, since due 
to the very high binding energies of monopolium, asymptotic monopoles might only be 
found, for short periods of time, in the center of galaxies, or clusters of galaxies, not in 
our laboratories. 

In the present work we aim at determining the existence and the dynamics of monopoles 
in the laboratory. The philosophy behind the present calculation is that the standard 
model allows for the existence of monopoles which are spin bosons. Therefore, given 
the appropriate kinematics we should be able to produce them. However, past experience 
has taught us that this is not feasible because, most probably, their mass is very large 
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and their production is outside our present experimental capabilities. Our proposal here, 
is that due to the large coupling between monopole and anti-monopole the two bind to 
form a low mass monopolium state. This state can be produce as an intermediate virtual 
state and we study its subsequent decays. Thus, in an indirect way monopole physics can 
be revealed. 



2 Monopolium detection 

We proceed to discuss signatures of monopolium, the monopole-antimonopole bound 
state, when produced in e + e~ annihilation^- We use to describe the interaction the low 
energy effective theory of Ginzburg and Schiller [23]. This theory is based on the stan- 
dard electroweak theory and in order to couple the monopoles to the photon and weak 
bosons one considers that m » mz , mz being the mass of the Z boson, and that 
the monopole interacts with the fundamental fields of the SU(2) ® U(l) theory before 
symmetry breaking, i.e., with the isoscalar field B, in the conventional notation of the 
standard model [25]. In this way the photon, 7, and the weak boson, Z , have the same 
coupling except for an additional tan9w , where 6w is the Weinberg angle, for the lat- 
ter. The effective description is based on the one loop approximation of the fundamental 
theory and therefore the effective coupling is proportional to g e ff ~ ^ g, where u is an 
energy scale which is below the monopole production threshold, thus rendering the theory 
perturbative. The dynamical scheme proposed by Ginzburg and Schiller leads to effective 
couplings in a vector like theory between the monopole and the photon [23], given by 

gj ff = C(J m ) g - = C(J m ) (2) 

with C(J m ) ~ 1, u the photon energy, m the monopole mass and N the monopole charge. 
The effective interaction between the monopole and the Zq becomes 

9eff = tan(0 w )0e// (3) 

where 6w is the Weinberg angle and naturally here oj refers to the Z energy. We have used 
the Dirac quantization condition Eq.(l) to express the coupling in terms of the electron 
charge. 

We study the process 



e + e- -> A^M + A' 

^B + C, (4) 

shown in Fig. 1, where A, A', B, C are 7's or Zo's, in all allowed combinations and M 
represents a monopolium state. We consider that the particle A' carries away the spin of 
the photon and therefore M represents the lowest scalar monopolium state. The coupling 
arises as a consequence of the generalization of scalar electrodynamics [23]. 

1 The description in terms of quark- antiquark annihilation is straightforward although complicated by 
the partonic description of the real experimental probes which are hadrons. 
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Figure 1: Diagrammatic description of the reactions studied. 



The standard expression for the cross section in these cases results in 
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Here M stands for the monopolium mass, P for the monopolium four-momentum and G 
is defined by 
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where Jm indicates the total angular momentum of monopolium. We recall that the so 
called unitarity bound, restricts the validity of the Ginzburg- Schiller approximation to 
M < m/Q [7]. The interesting physical situation occurs when M<m and consequently, 
from (121) and the fact that u ~ M, one gets g e ff ~ g(M/m) 1, which grants validity 
to the perturbative approach. 

We enter now the computation of the widths T. Taking into account that 
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where m A is the mass of the A particle in Fig. 1. Here q = (2 E e , 0). One obtains, 
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Note that the monopolium mass could be small and therefore we have to keep the mass 
term in the denominator. 

We now proceed to calculate 

T bc = 2tx\M bc \ 2 P {E BC ) 

Using the standard result [26] we obtain 
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where the approximation m ^> mB } mc has been used. Here m# and rric represent 
respectively the masses of the B and C particles of Fig. 1. 

The cross section, generalized to include the propagation of an unstable particle of 
width Ta (for us A = Z ), becomes 
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where Co indicates the number of Zq's present among particles A, A', B and C and we 
have used energy-momentum conservation at the vertex. 

In order to go ahead with the calculation, one has to obtain the wave function corre- 
sponding to monopolium. This is done and analyzed in the next sections. 



3 Monopolium Potential 

Our calculation of monopolium reduces to a quantum mechanical bound state calculation 
which provides us with its mass and its wave function in the relative frame. We will use a 
static non relativistic approximation and therefore our first step is to define the potential 
that binds the poles to form monopolium. 

We restrict our calculation to the lowest charge monopole, i.e. N = 1 in the Dirac 
condition Eq.(l). We regard the monopole as possessing some spatial extension in line with 
the arguments of Schiff and Goebel [271 [28] . This assumption makes the potential energy 
of the monopole-antimonopole interaction non-singular when the relative separation goes 
to zero. Mathematically we describe this feature by means of an exponential cut-off in 
the interaction potential, 

v[r) = - f p-«p(-/"-) j . (9) 

We fix the cut-off parameter /i by physical arguments. Eq.([H]) has the following prop- 
erties, 

i) r — > oo 

V(r) -> -— (10) 
r 

ii) r — » 

V(r) - + (11) 



When the monopole-antimonopole are closest to each other the distance between the 
corresponding centers O and O' is Too 1 ~ 2r m , where r m is the pole radius. For our 
estimates r m ~ 4:r c i assica i seems reasonable since this choice will allow the monopolium 
bound state to have very small mass for strong binding. Then 
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Figure 2: The figure on the left shows the Coulomb and the Coulomb with exponential 
cut-off potentials. The figure on the right shows energy levels of the Coulomb potential 
used as energy levels of the cut-off potential. Note that the lowest energy states of the 
cut-off potential correspond to excited states of the Coulomb potential. 
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T 

Consequently, the effective potential finally becomes 

1 - exp f-2 — r - — ) 

V(r) = ~g 2 " 'classical ) % Q3) 

r 

Note that with our choice, for r — > 0, V(r) — > —2m. Thus, the mass of the bound state 
becomes the energy over the minimum 

M = 2m + E binding . (14) 

Summarizing, our analysis shows that the cut-off potential is quite close to the Coulomb 
potential as long as the monopole radius, r m is greater than the classical monopole radius 
t 'classical- Thus, we shall use the "magnetic" Coulomb potential (Fig. 2) as our interac- 
tion in what follows. However, it is important to note, as will be shown next, that the 
lowest energy states of the cut-off potential correspond to excited states of the magnetic 
Coulomb potential (Fig. 2.). 

We use a non-relativistic approximation whose validity we will shortly discuss. Solving 
the Schrodinger equation for monopolium we obtain its binding energy, and therefore the 
mass of the system is given by [29] 

M = 2m-(^-) 2 ^>0. (15) 

where a = e 2 = 1/137 and n is the principal quantum number. We see that we can reach 
zero mass for n ~ 12 and therefore for n > 17 the formula is well defined and describes 
all values of M 
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Figure 3: Mass of the monopolium as a function of the size parameter. The unitarity 
bound corresponding to p ~ 0.41 is shown. 



< M <2m 



The monopolium radius is given by 

= 48aV. (16) 

^classical 

Now we introduce the size parameter p = r M /r dassica i. 

By substituting n 2 from Eq. fll6j) into Eq. (115|) . we obtain an equation for the monopo- 
lium mass as a function of its size, namely 

M = m( 2 -1), (17) 

which is plotted in Fig. 3. Although for low values of p our approximation becomes worse, 
we expect, that the soft behavior of the wave function at the origin, allows for order of 
magnitude estimates. This formula is extremely important in our development because 
it transmutates principal quantum numbers of the Coulomb potential into mass scales 
which are crucial to substantiate our scenario. 

The perturbative expansion of Ginzburg and Schiller used in the calculation of the 
production process limits, due to unitariry, the maximum value of the ratio of — < 1/6 
[TJ. We incorporate this external additional requirement to the bound state calculation, 
which is free from it, by showing the bound in all relevant figures. 

Before we continue, a discussion on the validity of the non-relativistic approximation is 
necessary. Let us therefore analyze how relativistic corrections will affect our calculation. 
As the monopoles are spin bosons one may describe the system by a Schrodinger type 
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Figure 4: Shown are the non relativistic 3 (solid), the relativistic bosonic 3 re i (dashed), 
the first order correction 3 re i (dot-dashed) and the ratio of an upper bound to the order 
/3 potential term to the binding energy to order (3 A with correct kinematics (dotted). The 
unitarity bound is also shown. 



equation and the first corrections to it do not start at order 3 2 but at order /5 4 (apart 
from merely the kinetic terms) [30]. These additional terms to order /5 6 are, 

A 4 A 6 g 2 rA2 rA2 1 

" -[A 2 , [A 2 ,-]]. (18) 



8m 4 16m 6 ' 32m 5 ' ' r 
They can be treated as perturbations to the Schrodinger equation. In the appendix we 
give detailed account of the calculation of these corrections. 

We show in Fig. 4 that the correction to the potential energy can be neglected safely, 
for an order of magnitude estimate, if the kinetic terms are fully taken into account. 
Moreover, we show that the minimal relativistic correction, consisting in the simplest 
approximation to the Klein-Gordon equation, namely in incorporating the mass [30] . 
leads to a velocity 

= TJ^, (19) 

1 ' 2 

which is almost indistinguishable from the more complete given by Eq. (l34l) . 

Furthermore, since our potential is cut off for small values of r we expect a slow down 
of the particles with respect to the conventional Coulomb potential and therefore the 
non-relativistic treatment is more accurate than in the pure Coulomb case. Moreover, the 
calculation for the wave function at the origin is less sensitive to the short range behavior 
of the potential, than the velocity which depends on the slope of the wave function. 
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4 Cross section estimates 



We proceed to use the formalism just developed to describe the production and decay 
of monopolium. The analysis that follows is physically appealing because the mass of 
monopolium may be chosen small, much smaller than the monopole mass, thus detection 
can occur at relatively low energies. Monopolium production is accompanied by radiation, 
which is also described by the formalism. Furthermore the calculation is easy to perform 
and physically understandable. 

It seems therefore safe to go ahead to calculate the monopolium decay probability as 
a function of p. The range of values of p : 3/8 < p < oo. Moreover, 

and therefore, given a value of p, one can determine n and this fixes |t/>(0)| 2 , which is what 
one needs for computing the decay probability. In summary, the calculation seems to be 
feasible in terms of only one mass scale, the mass of the monopole, m, and one parameter, 
P- 

Let us consider the case when monopolium is produced in its ground state, its wave 
function will have I = 0. Consequently [29] 
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one has 

1^,0,0(0)1 = ^-4= ^ 



The reduced mass of the monopolium system is m/2 and the Dirac condition Eq.(l) 
can be written as 

2 2 1 1 

9 e = - => a g a e = - 

one gets 

1 / m \ 3 / 2 
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Finally one can write the wave function in terms of the variable p to obtain 



... if 2y/3m 
|fc,o,o(0)| = - [ — 



3/2 



(23) 



This is the main ingredient to be included in the expression for the cross section that was 
computed before. 

For spin zero monopoles and being particularly interested in n large and i = , one 
has G ~ 1/4. Replacing the value of the wave function given in (1231) . and using the 
relation between M and p the cross section, Eq.flBJ), becomes 



a 



1 27 (tan6 w ) 2Nz o M 2 {4 f 2 _ M\ 



4 6 32 7T £ 



cr 



77T 



m J 



E\, E 2 B E 2 C (2 E e 



\{AE 2 e - m\Y + {m A Y\\ ((2E e - E A ,) 2 - E% - M 2 ) + M*T 



M 



(24) 



We proceed to describe the case when A is a photon propagator and the external 
particles A', B and C are outgoing photons. Clearly A = Z , A' = B = j,C = 7 also 
contributes to the 37 cross-section. We omit it here for simplicity, since we are just 
estimating the observability of the process, not its precise magnitude. 

In the case under consideration the above cross section becomes, 



4 6 7T 5 Q! 3 m 4 



M 2 / , M\ 3 



rp 2 7712 7712 

£j A , £j b £j c 



m J 4E 2 



((2E e - E A ,) 2 - E% - M 2 y + M 2 T 



M 



(25) 



It is evident from the above equation that the cross-section has a resonant structure for 



E A t — E e 1 



M 
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In order to maximize the cross-section we will set ourselves on top of the resonant peak 
and moreover will choose 



E B = E C = ^ [ 1 



M \ 2 \ 

2Ej r 



which are the values for the energy that maximize the numerator, E B Eg, while satisfying 
the conservation equations 



Eai + En + E C = 2E e . 



and 



PA' + Pb + Pc = 
These substitutions lead finally to the formula 
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The monopolium width is dominated by the 2 7 decay, 



Y M ~ T BC (E B = E C = M/2) 




(27) 



with a correction from the Zq decays which is about 30% and which for the purposes of 
our calculation is irrelevant. Using this width we obtain for the cross section, 



The results of our calculation are shown in Figs. 5 and 6 which we now comment. We 
have two parameters in our calculation, namely the monopole mass m and the monopolium 
mass M. We use in the plots M and the ratio M/m. Unitarity imposes a restriction on 
the latter M/m < 0.15 as can be seen in Fig. 3. We show data for two values of this 
ratio, 0.01, which serves for the purpose of developing our idea and satisfies clearly the 
unitarity bound, and 0.0001, which could realize certain cosmological monopole scenarios. 
The equations are sufficiently simple for any one to play with them. In Fig. 5 we show 
on the left the cross section at resonance as a function on the mass of the monopolium 
and for a beam energy just above the monopolium threshold, i.e. 2 E e = M + 0.001 GeV. 
This extreme case leads to low values for the cross section due to the vicinity of the 
threshold zero (see for comparison the dotted line which is away from the threshold), but 
is physically very appealing because the photon of the first vertex carries almost no energy 
and we have a representation of the scalar monopolium decay with two photons appearing 
back to back with an energy of M/2. In this case the photon of the first vertex simply is 
present to carry away the spin of the intermediate photon but (almost) no energy. We see 
that the value of the cross section increases as Mj m decreases. Thus initially monopolium 
made of heavier monopoles, for a fixed mass, would be easier to see, if it were not because 
as shown in the right figure, its width decreases dramatically. 

We include the limits of the Tevatron and LHC, which in our way of presenting the data 
are trivial. Note however, that both in the Tevatron and in LHC, the expected processes 
are the inverse of the studied here, namely the 2 7 (and other mentioned alternatives) of 
producing monopolium [231 EI] • 

Let us emphasize at this point looking at the data the interest of what we have named a 
two scale scenario. The two scales are the monopole mass and the monopolium mass. Our 
scenario becomes very interesting if the two scales are very different. It is the existence 
of an additional scale, which is represented in our scheme by the parameter M/m, which 
renders the present investigation exciting. If monopolium is a strongly bound system, it is 
its relatively low mass which limits the usefulness of the accelerators for monopole physics 
and not that of the monopole, conventionally assumed to be very large. If monopoles could 
bind to an almost zero mass bound state we could study monopole physics at relatively 
low energies. 

In Fig. 6 we show for a fixed monopolium mass how the peak cross section changes 
with the energy of the beams. We see that the threshold effect is extremely narrow and 
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Figure 5: Logarithmic plot of the cross section at resonance (in fb) (left) and the monopolium 
width (in GeV) (right) as a function of the monopolium mass (in GeV). The vertical lines 
represent the Tevatron (M ~ 265GeF)[S] and the LHC (M ~ 1100 GeV) bounds. In the left 
plot the full lines have been calculated for a beam energy slightly larger than the monopolium 
mass, 2E e = M + 0.001, to avoid the threshold zero and maximize monopolium production. The 
dotted line corresponds to a value of E e WOMeV above the threshold. The different curves are 
obtained for the shown values of M/m = 0.01 or 0.0001. 
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Figure 6: Logarithmic plot of the cross section at resonance (in fb) (left) as a function of beam 
energy E e for a fixed mass of the monopolium of 1000 GeV . The different curves are obtained 
for the shown values of M/m = 0.01 : 0.0001. 



11 



that the cross section jumps immediately several orders of magnitude. Thus the three 
photon detection seems much more favorable than the two photon one. The signal is clear: 
one photon recoiling against two others, whose dynamics describes a resonant structure. 
Once the threshold effect disappears the cross section is relatively flat with energy. 
Let us summarize our findings. In our two scale scenario, 

i) the studied cross section has a resonant peak (see Eq.(32)) at the monopolium mass 
M; 

ii) the order of magnitude of the cross section is "almost" beam energy independent, 
once we are away from the threshold, and is consistent with observability in present 
day machines [231 EH E21 [33] for monopolium masses of up to 1000 GeV (see Figs. 5 
and 6) and monopole masses only limited by the validity of our theoretical approach; 

iii) a similar analysis can be carried out for 7 Z Q and 2Z decays; 

iv) a similar analysis can be carried out for hadronic production, complicated by the 
inclusion of the sub-structure of the intervening hadrons. 

v) Our analysis can generalized to study the 2 7, 7 Z and 2 Z monopolium production 
process which should be the dominating mechanism for LHC 

5 Conclusions 

We have performed an investigation looking for hints of the so far not seen monopoles. 
Our working assumption is that monopoles appear strongly bound forming monopolium, 
a monopole-antimonopole bound state, due to their strong electromagnetic interaction. 

We develop a scenario in which monopolium is produced and desintegrates into 27, 
7 Z and 2Z 's. We detail the structure and magnitude of the first of this processes to 
determine observability. We develop a two energy scale scenario, whose 

i) low scale is governed by monopolium and we consider for quantitative purposes that 
it could be reachable by present day machines. 

ii) and whose high energy scale is governed by the monopole mass and arises through 
the structure of monopolium. 

m » E e . 

Under these circumstances we can estimate the cross section as a function of the 
monopolium mass. The monopole mass is determined by the value of the cross section 
and any mass is attainable, the limitations only arising from the approximations used in 
our model. 

Since at present we can not calculate the monopolium parameters, M and Tm, the 
experimental endeavor is not easy. There are however some features which might simplify 
the task, 
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i) the resonance peak of the monopolium can be found in four exit channels 37, 27 Zq 
and 7 2 Z 's and 3 Z 's; 

ii) monopolium can be produced in an excited state before it annihilates, thus the 
annihilation process will be accompanied by a Rydberg radiation spectrum; 

iii) the same processes can be studied hadronically, the only complication arising from 
the inclusion of the hadron sub-structure. 

The calculated values for the cross sections, corresponding to reasonable monopolium 
mass scenarios, render our calculation interesting and this line of research worth pursuing. 
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A Appendix: Relativistic corrections 

Let us define the relativistic factor (3 = v/c through the equation 



(3 2 =< nlmim s \ — T:\nl mim s > . (29) 
wr 

It can be easily calculated using the exact expectation value to give [29J, 

(3 = If. (30) 

V 4 p 

This result, which coincides with the semiclassical treatment and the use of Ehrenfest's 
theorems, namely equating the centrifugal and Coulomb forces 
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This corresponds to equating the absolute values of the kinetic and the binding energies, 

Kinetic Energy = \Binding Energy\, (31) 

which gives 

9 , -1 n 2 
p ( 1 \ m 



m v 8 a J ri 2 
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that gives rise, using Eq. (jl6j) . to 

p 2 _ 3 1 
m 2 4 p 

Thus, the non-relativistic calculation is only truly valid for p > 3/4. 

One can easily incorporate the kinetic terms in Eq. lfTB"]) to the lowest order approxi- 
mation leading to, 

m —- — + — + ... 32 

v 2 8 16 ; K J 

Performing the virial theorem calculation, 

3 2 8 4 13 6 

= m ( 1 + Y-Y + W (33) 

and therefore the relativistic velocity turns out to be 

Prel = £2 _ £4 £6 ■ (34) 

2 8 16 

The non-relativistic velocity f3 and the relativistic one (3 Te i are shown in Fig. 4. Finally 
the correction to the potential compared to the binding energy becomes, 

3^4 < [A 2 , [A 2 , i]] > g 2 p3 ^ 137/33 



\Ebinding\ 32 m 3 128 

The upper bound is also shown in Fig. 4 
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